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Abstract 

Modeling the interaction between a non-uniform magnetic field and a rotating conductive object 
allows study of the drag force which is used in applications such as eddy current braking and 
linear induction motors as well as the transition to a repulsive force that is the basis for magnetic 
levitation systems. Here, we study the interaction between a non-uniform field generated by a 
cylindrical magnet and a rotating conductive sphere. Each eddy current in the sphere generates 
a magnetic field which in turn generates another eddy current, eventually feeding back on itself. 
A two step mathematics process is developed to find a closed form solution in terms of only two 
eddy currents. However, the complete solution requires decomposition of the magnetic field into a 
summation of spherical harmonics, making it more suitable for a graduate level electromagnetism 
lecture or lab. Finally, the forces associated with these currents are calculated and then verified 


experimentally. 


I. INTRODUCTION 


If a magnet is moved slowly along the surface of a conducting sheet, eddy currents are 
generated in the sheet that create drag forces, opposite the direction of motion, between 
the magnet and the currents. These forces are the basis for eddy current braking,’ where 
moving magnets are slowed down as they pass by a conductor, and for linear induction 
motors,” where a moving magnetic field can strongly propel a conducting sheet (see for 
example the electromagnetic aircraft launch system, EMALS).* Second order eddy currents 
are also induced that create an image magnet whose orientation opposes the real magnet. 
This repulsive force is insignificant at low speeds, but becomes dominant as the magnet is 
moved rapidly across the sheet, pushing the magnet away from the conducting sheet, i.e., 
providing increasing lift as the drag force diminishes. The most well-known examples of this 
are certain embodiments of magnetic railways, i.e., MagLev,* where a moving car floats on 


a conducting railway. 


This transition from drag to lift has been studied®’ but the derivations are typically ap- 
plied to infinite conducting sheets where an eddy current wake is created. These derivations 
are pertinent for many real-world applications, but they are concise, often only approximate, 
and complicated. In the discussion below, a step-by-step derivation will be provided showing 
how to calculate the drag and repulsion forces between a cylindrical magnet and a rotating 
solid sphere. We chose this geometry because it allows closed form solutions for all of the 
electro-magnetic entities involved and can yield insight into the physics behind this force 
evolution. Some of the details in this derivation are adapted from Heinrich Hertz’s doctoral 
thesis® where he develops a methodology for finding the eddy currents in a rotating sphere 
when exposed to a constant magnetic field. However, Hertz math is customized, and con- 
sequently difficult to follow, since the modern Heaviside notation had not yet occurred. So, 
from a historical perspective the development given below provides insight into an approach 


first developed by Hertz in the 1880s when he was only 23. 


The derivation provided in this paper is appropriate for a graduate level electromagnetism 
lecture or lab. The problem requires knowledge of many of the key concepts of electromag- 
netism and applies them with the goal of modeling and understanding important, but often 
neglected, topics such as eddy current generation and the interaction between a magnet and 


a conductor moving at some relative velocity. Our prior papers described the interaction of 


a magnet with a slowly rotating shell®, where torques were generated from the drag force 
but since only the first eddy currents were modeled, lift did not occur; and the interaction 
of a uniform magnetic field with a rapidly rotating solid sphere!®, where the symmetry of 
the uniform magnetic field prevents a lift force from occurring. This paper, describing the 
interaction of a non-uniform magnetic field, such as that generated by a cylindrical magnet 
or current loop, with a solid sphere rotating at any speed, does yield lift forces. Also, this 
is the simplest geometric form where a complete representation of the eddy currents can 
be obtained and where the results can be expressed in closed form. It should be noted 
that this complete representation requires the applied magnetic field be decomposed into a 
summation of spherical harmonic components. This two-dimensional orthogonal series rep- 
resentation may be beyond an undergraduate level, but offers an introduction into complete 
infinite bases, i.e., Hilbert Spaces, for the graduate student outside of those presented in 
quantum mechanics. In addition, this paper presents an experiment that can be performed 


to provide verification of the theoretical predictions. 


Il. STARTING CONDITIONS 


A sketch of the problem is shown in Figure 1 with both Cartesian and spherical coordi- 
nates that have their origin on the sphere center and are fixed to the laboratory frame of 
reference. The sphere rotates with a constant angular velocity given by © = w2 and has a 
radius given by R. It has uniform conductivity, 0, but has no net charge and is not magnetic, 
i.e., the relative permeability of the sphere is equal to one. Also, since the sphere is round, 
the conductivity, as seen from the lab frame, is not time varying (this is important since 
most rotating objects would cause the conductivity at some location in space to alternate as 
material entered and then left it). The applied magnetic field, B, is not time varying and is 
cylindrically symmetric since we are assuming a round magnet or current loop as the source 
of the field. The coordinate system has been rotated so that the magnetic field points in 
the x-direction and is symmetric about this axis. Since the magnetic field, angular velocity, 
and conductivity are not time varying, all deduced parameters such as currents, potential, 
and electric fields are also not time varying. 

We are only concerned with the magnetic field in the volume of the spinning sphere and 


assume that this applied field has been created by currents (or magnets) outside of the sphere. 
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FIG. 1. A magnetic field applied to a rotating solid sphere showing both the Cartesian and spherical 


coordinate systems used. 


Combining this assumption with the time independence of any electric fields causes Ampere’s 
Law to reduce to Vx B = 0. This implies that the magnetic field can be written as the 
gradient of a magnetic scalar potential,’ ®y,, ie., B=-V®vy. This potential is commonly 
used in magnetostatics because when combined with Maxwell’s equation, V - B= 0, it solves 
Laplace’s Equation, i.e., V?®,, = 0, a well-studied equation. In spherical coordinates this 
means the magnetic scalar potential can be written as a summation of two orthogonal sets 


of functions, Legendre functions and trigonometric functions,'! 


®x(r,0,¢) = yo S- Bam Fey PM (C08 0) cos(m@). (1) 


n=0 m=0 


This is not the most general form of the solution; polynomial terms that are unbounded at 
the origin have been dropped and, since the magnetic field is symmetric about the x-axis, 
only cosine terms are retained. Equation 1 is sometimes expressed in terms of Spherical 
Harmonics, but we chose to use the more explicit form shown. The factors in Eq. 1 have 
been chosen to simplify the analysis and to give the b,, coefficients the same units as the 
magnetic field, i.e., Tesla. Using this decomposition, the magnetic field can now be found 


allowing an analysis of the problem to begin. 
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Ill. THE EDDY CURRENTS 


It should be noted that while discussing eddy currents, we are really talking about current 
densities, but for simplicity will continue to refer to them as currents for the remainder of 
the paper. In this section, the primary task is to find the currents in the rotating sphere 
when exposed to an applied magnetic field. Start by using B = —V®, and the expression 
in Eq. 1, to find a representation for the magnetic field 

Be Batt (2) 
n=0 m=0 


where, 
—bnm (5) ae nP" (cos 0) cos(m@)r 


Bee +bnm ie ((1+n) cos @P™ (cos 0) + (m—1—n) P™., (cos 8)) csc 6 cos (m¢) 6 
+Bn mm (2)ns P(cos 0) csc 6 sin(md)¢d 


and ene is the n,m component of the magnetic field. A typical magnetic field of this form 


is shown in Fig. 2. Recall from Ohm’s Law that the current is related to the applied field 
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FIG. 2. A representative magnetic field created by a coil, aligned with the z-axis near the rotating 


sphere. 


by 
Tee (Eng +0x Be) (3) 


Using the angular velocity, W@ = w2Z, the velocity of any point in the sphere is given by 
C= x72 lo find ax Baym the magnetic field is converted to Cartesian form, the cross 


product taken, and the result converted back to spherical coordinates yielding 


0X Bum =— bash ((1+n) cos OP (cos 6) + (m—1—n)P™, (cos @)) cos(m@)F 


n (4) 


— bam wean sin OP" (cos 0) cos(m@)6 


Now, recall that conservation of charge states that the divergence of the current is equal 
to the change in the charge density with time, i.e., V- J= Op/0t. However, since all of the 
parameters in this problem are time independent this becomes V - J =0. So, taking the 


divergence of Eq. 3 and using Eq. 4 yields 


V+ Enm=—V- (0X Bam) 
n-1 (5) 


a ((1+2n) cos 0P"(cos @) + (m—1—n)P™., (cos 4)) cos(m@) 


= Damp 


The electric potential, V, is related to the electric field by E = —-VV so Eq. 5 becomes a 
Poisson Equation 
i Ane (@ x Brig) (6) 
If we can find the electric potential, then we can find the electric field and then the 
induced currents in the sphere created by the applied field. To do this we need boundary 
conditions on V and these are that the potential must be bounded throughout space, it must 
be continuous across the sphere outer boundary, and no currents can leave the sphere. This 
current boundary condition is physically reasonable—current cannot flow out of the sphere— 
and it is commonly used to solve eddy current problems. However, it is not immediately 
obvious that it is pertinent since the current is not explicit in Eq. 6. But, if we insist that 


the radial component of the current at the sphere boundary, J, (2) = 0, then Eq. 3 yields 
OV = 
(@ x B) (R). (7) 


o (E,(R) < (ax B) (R)) ates oat 
This results in a Neumann boundary condition on the potential (we have dropped the m,n 


Or 


r 


subscript to show the vector component index). Integrating the radial portion of Eq. 4, the 


electric potential inside of the sphere becomes 


wyrrtl 


Venn = Unt 
™ (n+1) Rr} 


((1+n) cos OP,” (cos @) + (m—1—n) Pl, (cos @)) cos(md) (8) 


and the homogenous solution of Eq. 6 is used outside of the sphere to ensure continuity of 
the potential across the sphere boundary. Using this form for the potential, the electric field 
can be found and, from Eq. 3, the current can be obtained; 


= rr 


Jnjm = — owb —"  escd (mrs (co 0) cos(md)O 


mm Rr-1 (n+1) (0) 


+ ((14n) cos 6P(cos 0) + (m—1—n)P”,, (cos 6)) sin(m}6) 


As first pointed out by Hertz, there are no radial current components; all of the current 
flows along shells, forming closed loops since the divergence is zero. Figure 3 shows a 
representative form for this current, looking along the x-axis. The velocity of the spheres 
surface shown in this Figure is in the positive y-direction and the magnetic field is into the 


paper yielding an upward force that results in the current shown. 
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FIG. 3. A typical form for the eddy current generated by the rotating sphere in the presence of 
a magnetic field aligned with the z-axis. This view is looking from the magnet in the positive 


z-direction with the origin at the sphere center. 


The eddy current shown in Eq. 9 is generated by the applied magnetic field interacting 
with the rotating conducting sphere, but it is not the total current. This eddy current 
generates a magnetic field within the sphere that, through the same process shown above, 


yields a second eddy current. Then the second eddy current generates another magnetic 


t 


field which yields a third eddy current, and, as will be shown below, this third eddy current 
is of the same form as the first eddy current but of opposite sign. So the third eddy current 
acts to cancel the eddy current in Eq. 9, yielding a negative feedback loop. 

To obtain the total current within the sphere, we will take the form of the current shown 
in Eq. 9 and generalize it by making the radial dependence an unknown function and the 
goal of the analysis will be to find this function. Once it is obtained, then the second eddy 
current can be found and the forces and torques on the sphere can be calculated. This 
insight, that the current shown in Eq. 9 can be generalized in this fashion—only changing 
the radial dependence—is due to the 23-year-old Hertz.® 

Let’s start with a general form of the first eddy current, de 


nym) 


caused by the applied field, Tees defined by Eq. 9, added to the third eddy current, Te 


which is equal to the current 


Fist da 


nym ~— 


= — 0W0nmCnm (rT) 7 esc 8 (mp (cos 0) cos (m@) 6 (10) 
+ ((1+n) cos @P"” (cos #) + (m—1—n) P,, (cos @)) sin (m@) 6). 


Cnm (1) is a dimensionless function of the radial distance and the goal of the analysis is 
to find this function. Note that c,.(r) and the other current related scalar functions 
derived below are functions of the subscripts n and m. If there is no third eddy current, 
then nm (r) = (m/(n+1)) (r/R)"™, which corresponds to the limiting solution when the 
sphere is small or rotates slowly and must emerge from the more general solution. The 
process will be to find the magnetic field generated by the current Fao and then find the 


second order eddy current, Je generated by that field. We then do this again to find the 


third order eddy current, which will allow the feedback loop to be closed and the function 
Cnym (r) to be found. 


In order to find the magnetic field Be, generated by the current, J‘ use the vector 


nym) 


potential equation, VAL = p19 F mas (where ju9 is the magnetic constant, 4a x 1077 A/m/?) 
to find A? ,, and then use B? ,, = V x A?.,, to find the field. Note that we use superscript “2” 
here to indicate that these are the second set of magnetic fields and potentials and not that 
these entities are squared. To solve the vector potential equation, the current and the vector 
potential are first both put into the form J(r, 0, ¢) = Jy (7,0, 0) &+ Jy (7,0, &) G+ Jz (7, 4, d) Z 
so that the Laplacian operator can be applied to each component of the vector functions, 


but the dependence is still in spherical coordinates so that the boundary conditions can be 
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easily handled. Doing this to the current yields 


Phy = Fuh nn Cam (P) (5 ((n+n?-+m—m?) Pe (cos8) cos ((m=1) 8) 


+ Pi" (cos @) cos ((m+1) 6))# 


| 5(- (n-+n?+m—m’) P (cos @) sin ((m—1) 6) + P?”** (cos 6) sin ((m+1) ¢))9 


+mP"" (cos 6) cos(md) 2). (11) 


This expression shows that in a Cartesian decomposition the current components are all in 
the forms of spherical harmonics (note the order of the trigonometric functions matches the 
superscript of the associated Legendre functions). This is important because the Laplacian 
of a spherical harmonic returns a modified radial function times the same spherical harmonic. 
So, the inhomogeneous portion of the vector scalar potential can be put into the same form 


as the current but with a different radial dependence, i.e., 


At = Cra dager (5 ((n+n?+m—m’)P”~ "(cos 8) cos((m—1)¢) 


+ P™** (cos) cos((m+1)¢))4 


; (—(n tn? +m—m?) PP" (cos 6) sin((m—1)¢) + P™*' (cos @) sin((m+1)¢)) g 


+ mP""(cos 0) cos(m¢)3), (12) 


where the new dimensionless radial function, gpm (7), is given by 


roan) + Ar gi m(T) + (2—n—n’) Inm (r) +r? pgowe(r) = 0. (13) 


However, the form for the vector potential shown in Eq. 12 is not complete. The homoge- 


neous solutions must be added in order to obtain continuity of each component of the vector 


=> 


and each component of the magnetic field, B? 


nym) 


potential, Aes across the sphere boundary 
where r = R. Since each spherical harmonic in Eq. 12 yields a corresponding homogeneous 
spherical harmonic both inside and outside of the sphere, there are ten coefficients to be 
found and determining them is a tedious process. In addition, Eq. 13 for gn (7) has two 
homogenous solutions, one of which is determined by stating that gn m (7) must be bounded 
at the origin. The other, when placed into the complete solution for AP in drops out, having 
no effect on A2_, so we are free to choose an arbitrary form for this second homogeneous 


nym? 


solution which simplifies the resulting form for both AD and for gnm (r). The final result 
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for the vector potential in spherical coordinates is 


Drm Inm(r) rmesc OP™(cos 4) cos(me) 4 


( 

ae Drm Inm(r) r (mcot OP?"(cos 0) + P™*!(cos 8)) sin(m@) 6, PeOR 
a brish Gum Ry Re msc 6P™(cos @) cos(me) 6 

( 


+Bnm Gnjm(R) = = (mcot P™(cos 0) + P™*(cos 0)) sin(md)¢, r>R. 
(14) 
and for gnm(r) is 
f capt 1 enm( ti 
Gnm(7) = pow |r? / a Sneae soree i, a+ mn) dt | . (15) 
r 0 


We have converted the vector potential back to a spherical basis so that curl can be taken 


to yield the second magnetic field, Be es For the volume within the sphere this is given by 


Be = —Dpm M142) gn m(r)P2"(cos 8) sin(m®) F 
+ bam (29nm(?) + 79hm(r)) esc @ ((1+n) cos 6P."(cos 8) + (m—1—n) P™,; (cos 4) sin(md) 6 


— Drm ™ (29nm(1) +r In (P)) csc OP" (cos 8) cos(m@) db (16) 


Comparing this to the applied field, the radial dependence is of course different, but the 
angular dependence is identical except that cos(m@) has become — sin(m@) and sin(m@) 
has become cos(m@) indicating a net rotation of the field about the z-axis. 

Following the same procedure shown above for the applied magnetic field, we next need 
to find v x Be: So, convert the magnetic field in Eq. 16 to Cartesian form, take the cross 


product, and convert the result back to spherical coordinates yielding 


UX Bet. = —DpmWt (2Gnm(") + 19nm (1) ((n+1) cos 6P"™ (cos 6) 
+ (m—1—n) P™., (cos @)) sin(md)F 


n 


— nm Inm(r)n (1+n) sin 6P”(cos 6) sin(m@)é. (17) 


The potential inside of the sphere is obtained by integrating the radial portion of Eq. 17, 


using the current boundary condition, yielding 
ges = —bn mT Gnm(T) ((1+n) cos 0P)" (cos @) + (m—1—n) P77, (cos 0)) sin(m@), (18) 
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which is very similar to the potential generated by the applied field shown in Eq 8. 
The electric field can be found from the potential and, combining this with the cross 


product result from Eq. 17, the second eddy current (see Ohm’s Law in Eq. 3) is obtained: 


ioe = —owm nm! Gnm(") csc OP™ (cos 9) sin(m@) 6 


+ awmbnm Inm(1r) esc (1+) cos 6P;" (cos 8) + (m—1—n)P7", (cos #)) cos(md) ob. (19) 


The impressive result here is that the second eddy current is nearly identical to the first 
current shown in Eq. 10. The function c¢,,(r) has become Mgr m(r) and the ¢ dependence 
has shifted; cos(m@) has become sin(m@) and sin(m@) has become —cos(m@), indicating 
a rotation of the vector field about the z-axis. Figure 4 shows a representative second 
eddy current. The circular shape taken by this second eddy current causes it to generate a 


magnetic field that opposes the applied magnetic field. 
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FIG. 4. A typical form for the second eddy current generated by the rotating sphere in the presence 
of a magnetic field aligned with the z-axis. This view is looking from the magnet in the positive 


z-direction with the origin at the sphere center. 


In order to find the third eddy current, Te iy start with the current shown in Eq. 19 
and repeat the process from Eq. 10 through 19. The mgp.m(7) will become m?ky m(r), the 
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cos(m@) will become sin(m@), and the sin(m@) will become — cos(m@), yielding 


Aker = Twby mM knm(r)r csc 0P’” (cos @) cos(m@) 6 
+ aWbym mM kn m(r)r esc 6 ((1+n) cos OP’"(cos @) + (m—1—n)Pi", (cos 0)) sin(m@) b (20) 


and 
R r 


ee Inm{t) 1 / Inm(t)t"* 
= | : 21 
Kenn (r) HMogwW | r / (1 al Wnjir-? dt pnt2 (1 + 2n) dt ( ) 


r 


This third eddy current has nearly the same form as the first eddy current (Eq. 10); the only 
significant difference being a sign change, indicating negative feedback, which is needed to 
achieve a stable solution. 

Recall that the starting current is the sum of the current created by the applied field and 
the third eddy current, i.e., inane = bane + Teas So, combining Eqs. 9, 10, and 20 we obtain 


Cnm(?) = (m/(n+1))(r/R)" — mK nmr). (22) 


The functions Chm(7), Gnm(7), and Knm(r) are also related by the second order differential 


equations (see Eq. 13) 


1Gnm' (1) + Ar Gam (7) + (2—n=N")Gnm(r) +r? Moowe(r) =0 

Tham’ (r) + 4rknm'(r) + (2-n=n7)knm(r) + 7? upowg(r) =0. (23) 
Solving Eq 22 for the function k,,,,(r) and substituting this into the second equation in 
Eq. 23 yields 

Cam (tT) + Arn (r)+(2—n=—17 )enm (7) —1? gowm? gn m(r) =0. (24) 
Now Eq. 24 can be solved for gn(7) and substituted into the first expression in Eq. 23 to 
yield a fourth order differential equation in (17) 
r'Cnm (r) +8r36n m?) (r)—2(-6+ nF n)r? en mi" (r)—4n(1+n)renm' (r) 


+(-2n—n?+2n?+n*4 +am'r*)enm(r) =0, (25) 


where a?, = figowm and a» has units of inverse length. This fourth order differential 
equation has four solutions corresponding to each of the Kelvin functions, but two of these 


are unbounded at the origin and will be dropped since we require bounded solutions at the 
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origin. However, for the case of a thick walled rotating shell with no material at the origin all 
four solutions would be required to meet the boundary conditions. So, the general solution 


to Eq. 25 is 
Cra Tt) = (Cini belts Om?) +Conmbeing1/2(Omr))/(Amr)?”. (26) 


Using the second equation in Eq. 24 g(r) can be found 


g(r) = (1/m)(Canmbetn41/2(Am?) — C1,nmbein41/2 (Aim?))/(Omr)?. (27) 


Since the general solutions shown in Eq. 26 and Eq. 27 simply yield each other when sub- 
stituted into the differential equations, the coefficients cannot be found using the differential 
relationships. Instead, the integral relations, Eq. 15 and Eq. 21, along with Eq. 22, are used. 
The expressions shown in Eq 26 and Eq 27 are substituted into these integral relations and 
two linear expressions for the coefficients, cjand cz, are obtained. Solving these yields the 
following expressions for the coefficients, 


(H2n)Vam R(bein 1/2 (am R)+ber,_1 /2 (Om R)) 


C1nzm = — V2(re+1) (bei? j9(mR)+ber?, 4 /9(amR)) 
é es (142n) Vem R( bein /2(@m R)—bern1/2(amR)) (28) 
arm V2(nH1) (bei, (am R)+ber2_, 79(amR)) 


Combining Eqs. 26-28 and recalling that a, = ,/figqwm yields the two eddy currents, ye 


nym 
72 
and J 


nym) 


from Eqs. 10 and 19. The sum of these two current densities over all values of n 


and m yields the total current density, tine in the sphere, where 


n=0 m=0 
T= Si 
n=0 m=0 
fe= JP, (29) 


IV. THE NET FORCE 


The primary goal of this paper is to obtain the drag and lift forces between the source coil 
or magnet and the rotating sphere. However, there are several sources of force distributions 
within the sphere and only a subset yields a non-zero net force. For example, as shown 
above, the rotation of the sphere in the presence of the applied magnetic field yields forces 


that cause currents to flow. However, since there is no net charge on the sphere, for every 
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force on an electron there is an equal and opposite force on a proton. So, these forces exert 
no net force on the sphere. 

Another source of force is the interaction between the currents themselves. Every dif- 
ferential current element in the sphere feels a force from its interaction with every other 
differential current element. However, between any two current elements this force is equal 
and opposite, i.e., the force one current element generates on a second current element is 
equal to, but opposite, the force the second current element generates on the first one. So 
if the current elements are both inside of the sphere their net force on the sphere is zero. 
This is not to say that these forces are small. The internal stress inside of a current can be 
equated to those inside of a pressure vessel. But the net force between currents J* and J: ei 
or between themselves, is zero because they are all located within the sphere. 

The only remaining source of force is the interaction between the current that generated 
the applied field and the currents inside of the sphere. As before, this interaction yields 
equal and opposite forces, but one acts on the sphere and one on the coil or magnet that 
generated the applied field, so these forces affect separate entities that can respond to that 
force. In practice, this interaction is usually described through the magnetic field generated 
by one of the currents and the choice is arbitrary. In the present case, the currents and fields 
inside of the sphere are known and the sphere has useful symmetry so we choose to use the 
applied field, B , inside of the sphere and the currents, J and J”. So, the total net force on 
the sphere, with an equal but opposite force on the magnet is given by the Lorentz Force 
Law: 

pe i i Bae (30) 
sphere 

The total force can be calculated by finding the summations in Eqs. 2, 10, 19, and 29 and 
substituting the results into Eq. 30. It is recommended that the reader convert each of the 
vectors over to a Cartesian basis but maintain spherical coordinates as in Eq. 11 to simplify 
the computations. With the exception of the radial integral, much of the force calculation 
can be done symbolically, but it is a long tedious process with little physical insight to be 
gained. Our recommended approach is to decompose the incident magnetic field to find the 
coefficients, by.m (as discussed in the next section) and substitute them into the expressions 
to find the forces numerically. The first current, J‘, will interact with the incoming field 


and produce a force component that is only in the y-direction. So, if only the first current 
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is present then the magnet will feel a force pulling it along in the direction of the spheres 
velocity, while the sphere will feel an equal but opposite force. The force on the sphere will 
be concentrated near the magnet and, since it is not directed through the center of mass of 
the sphere, will generate a torque that will try to slow down the sphere. The force due to 
the second eddy current, BP, is directed only along the x-direction. It acts as a repulsive 


force, pushing the magnet away from the rotating sphere. 


V. MAGNETIC FIELD DECOMPOSITION 


In order to calculate the forces, the coefficients b,, must be found. We chose to use the 
radial portion of the magnetic field expression shown in Eq. 2, rather than the magnetic 


scalar potential, because the magnetic field is often known and is continuous. 


=e 3 3 Dia (Z)" npr (cos 0) cos(m¢@) (31) 


n=0 m=0 
Start by multiplying both sides of Eq. 31 by cos(m’¢) where m’ is an integer between 0 and 
n and integrate ¢ from 0 to 27. This integral is zero unless m = m’, in which case it is equal 
to either 27 when m = 0 or 7 when m 4 0. 
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/ SCR GO Hm) Youn (= -) nre(cos8), (32) 


where v(m) = 1 ifm ¥ 0 and equals 2 if m = 0. Now multiply both sides by P”(cos 6) sin 6 
where n’ is a positive integer and integrate 6 from 0 to 7. The integral on the right side is 


zero unless n = n’ in which case 


| [om mo)B,(r, 6, ¢) P’"(cos 0) sin 6 dé dé 


r ‘- 2n(n+m)! (33) 


~ Han) gene 


T2 


Equation 33 is true for any value of r so we choose r = rg and then rearrange to solve for 


the coefficients, yielding 


be See m)! 


sf r,0,@) P’*’ (cos @) sin @ cos(m@) d@ dé. (34) 
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Note that the b,, coefficients must be finite son 4 0. Also, any bp» that appear in the 
equations above that have negative n or m values or m great than n are defined to be equal 


to zero. 


VI. THE EXPERIMENT 


In order to test the theory, we positioned a cylindrical rare-earth magnet near a rotating 
solid sphere and measured the two force components. The sphere was 0.051 m (2 inches) in 
radius and composed of aluminum with a conductivity of 2.3 x 107 S/m. It was mounted 
with bearings on both the top and bottom and rotated up to 60 Hz. We strongly encourage 
eye protection be worn when the sphere is rotating and recommend a protective shield or 
cage be placed over the sphere and motor. Slight misalignments can cause resonances that 
may lead to bearing failure. 

A black cloth was attached to the sphere, so that a light source and detector could monitor 
the reflectivity of the sphere and determine its rotational speed. A force gage (Imada DS2-1) 
was used to measure the lateral, or drag, force and the repulsive, or lift, force (see Figure 
2 for the configuration used to measure the repulsive force). The rare earth magnet was 
0.025 m (1 inch) in radius and 0.013 m (0.5 inches) thick and the center of its face was 
located 0.02 m (0.79 inches) from the spheres surface. The magnet’s field was measured 
using a Gauss probe and since the magnet is relatively thin, it was modeled as a round 
current loop!® with an equivalent current of 9600 Amps. The magnetic field of a round 
current loop is given in Eq. A4, but in order to use this equation it had to be converted to 
Cartesian form, rotated so the loop would be aligned with the x-axis, shifted back away from 
the sphere, and then converted back to spherical coordinates so that the b,,, components 
in Eq. 34 could be found. 

Data was taken and compared with the theoretical model discussed above and the results 
are provided in Fig. 6 and Fig. 7. The match is remarkably good below 30 Hz but begins to 
deviate above this frequency. One explanation for the variation is that both forces are present 
on the sensor, even though we are only measuring one axis at a time. While measuring the 
repulsive force, a side load exists that peaks in the 30-60 Hz range causing some flex in the 
shaft and changing the alignment of the magnet slightly. A similar effect is present in the 


measurement of the lateral force, where the repulsive force continues to grow and will push 
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FIG. 5. The experimental configuration used to measure the repulsive force between a spinning 


solid sphere and a rare earth magnet. 


the magnet further away from the sphere, thereby decreasing the lateral force measured 
by the sensor. Additionally, the sphere experiences a resonance when spinning up through 
30-50 Hz that causes quite a bit of vibration and movement of the sphere. This is most 


noticeable in the data set for the lateral force measurement. 


The models were fit by allowing the distance between the magnet and the sphere to vary. 
The best fit for the repulsive force model occurred when the magnet was assumed to be 
0.017 m away versus the 0.02 m as configured. Likewise, the lateral force model gave a best 
fit for a distance of 0.021 m versus the 0.02 m as configured. A straight edged cardboard 
shim was used to adjust the distance between the sphere and magnet for the test. Its not 
unreasonable that there may be a few millimeters of error in the alignment, since perfect 
alignment of the center of the magnet with the outer edge of a curved surface is difficult to 
achieve. Additionally, the side loads and vibrations also vary the alignment during the test. 
Another potential source of error is that close up, the neodymium magnets magnetic field 


will vary some from the idealized coil model we used to represent the magnet. 


Overall, the data closely agrees with theory. However, Figures 6 and 7 don’t show the 
full extent of the variations in the forces with increased rotational frequency. In fact, as 
discussed in our previous paper, the forces are actually dependent on a combination of the 
sphere’s radius, electrical conductivity, magnetic permeability, and the rotational frequency. 


The theoretical forces for our test sphere are plotted below in Figure 8 as a function of 
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FIG. 6. Data (blue) and theory (black) showing the repulsive force. 
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FIG. 7. Data (blue) and theory (black) showing the lateral force. 


the dimensionless parameter, r2,/ji90w, to better illustrate the full range of force behavior. 
The first eddy current expression is responsible for generating the lateral drag force on the 


magnet while the sphere experiences an equal but opposite force in the y-direction. This 
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force climbs and peaks and then drops off again with an increase in r2,/fioow. Additionally, 
the non-uniformity in the force density across the sphere results in a decelerating torque. The 
second eddy current generates the repulsive force between the magnet and the sphere and 
continues to climb with r2,/pioow, but levels off as the incoming magnetic field is becomes 


increasingly rejected inside the sphere by the described negative feedback process. 
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FIG. 8. These plots show the changes in the lateral and repulsive forces with an increase in 
ro,/p97W. For these plots, r2, Wo, and o were assumed to be constants associated with our test 


sphere. 


VII. CONCLUSION 


It has been shown that the eddy currents generated by a fixed cylindrical magnet oriented 
along the equator of a spinning solid conducting sphere can be represented by two closed form 
expressions. These expressions where found by solving a fourth order differential equation 
with appropriate boundary conditions, yielding a solution in terms of Kelvin functions. The 
experiment showed close agreement between the measured forces and those computed with 
the derived expressions, confirming our revised version of the theory first worked out by 


Hertz. 
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Appendix A: THE MAGNETIC SCALAR POTENTIAL 


The magnetic scalar potential, ®),, arises in regions of space where there are no time vary- 
ing electric fields or source currents, allowing the magnetic field to be written as B=-V6,y, 
where V2®,, = 0. As with the electric scalar potential, manipulating scalar functions is 
easier than vector valued functions, so one might suppose that the magnetic scalar potential 
would be commonly used. However, it is somewhat difficult to calculate for even a simple 
current loop and, as pointed out below, it is not a continuous function. 

It can be shown that the magnetic scalar potential is given by the current in a closed 


loop times the solid angle of that loop as seen from a location in space. Analytically 


m= pol f Lae (Al) 
A 


where J is the loop current, A is the area of the loop, 7 is a unit normal out of the loop, 
and d is a vector connecting the observation point to a location on the surface of the loop,!° 
This interesting geometrical interpretation was known to Maxwell, who provided a form 
for the magnetic scalar potential for a round current loop, see Fig. 9, in terms of a series 
expansion! and after him other series expansions for the current loop were generated.!:! 
At that time, it was preferable to represent functions as series expansions that could be 
calculated using look-up tables as opposed to using functions such as elliptical integrals. 
However, with modern computers closed form functional representations can be handled 
easily and are preferable. A more recent publication!” provides such a solution in cylindrical 
coordinates. It is presented below in spherical coordinates. 


Define two parameters 


A4ar sin 6 Aar sin 0 


p= q) fe 
@+r+2arsind © p (a +rsin0) 


(A2) 


where a is the radius of a current loop. The magnetic scalar potential for a round current 


loop, with current J, in spherical coordinates is 


b pot /1 —k/6( K(k) — V1 — GII(6,k)) when a<rsiné (A3) 
Hol |/1 — k/B(K(k) + /T— BIL(8, k)) when a>rsin@ 
where K(k) and II(6,k) are complete elliptical integrals of the first and third kind, re- 


spectively. It should be pointed out that the scalar magnetic potential is not a continuous 
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viewpoint 


FIG. 9. A current loop showing a viewpoint location. 


function. The magnetic field follows closed paths since the field lines must end on themselves, 
indicating that the magnetic scalar potential has positive (or negative) slope at all locations 
along one of these paths. However, a function cannot continuously increase (or decrease) 
along a closed path and end up with the same value. The resolution to this quandary is to 
first note that the magnetic scalar potential form shown in Eq. A3 is missing an arbitrary 
constant and that the sign (due to the square roots) is ambiguous. By adding constants 
and changing the sign of the function it can be made continuous, steadily increasing (or 
decreasing) until one passes through the area of the loop. Across this area, in the x/y plane, 
the solid angle of the viewer changes from 27 to —2z (or vice versa) reflecting the change 
in the sign of the dot product in Eq. Al. This discontinuity allows the magnetic field to 


continuously follow a closed path. 


It can be verified with programs, such as Mathematica, that the gradient of Eqn. A-3, 


2 


yields, to within a sign, the magnetic field for a current loop [18]; 


~ acos@ . 
B= Lol yr 
n(a2+r?—2ar sin 0)Va2+r?4 2ar sin 0 
csc 0 “ 
+ Mol r?+a? cos(20)) E(k) —(a?+r?—2ar sin 0) K(k)) 0 
ne 2n(a2+r?—2ar sin 0)V/a2+r2+4 2ar sin 0 (( EES YU) 


(A4) 


where E(k) is the complete elliptical integral of the second kind. 


22 


10 


11 


12 


13 


Mark.A.Nurge@nasa.gov, (corresponding author) 

Robert.C. Yougquist @Qnasa.gov 

Wikipedia, “Eddy current brake — Wikipedia, The Free Encyclopedia,” 2017, https://en. 
wikipedia. org/w/index.php?title=Eddy_current_brake&oldid=772624886. 

Wikipedia, “Linear induction motor — Wikipedia, The Free Encyclopedia,” 2017, https:// 
en.wikipedia. org/w/index.php?title=Linear_induction_motor&oldid=775957988. 
Wikipedia, “Electromagnetic Aircraft Launch System — Wikipedia, The Free Encyclope- 
dia,” 2017, https://en.wikipedia.org/w/index.php?title=Electromagnetic_Aircraft_ 
Launch_Systemkoldid=769672993. 

Wikipedia, “Maglev — Wikipedia, The Free Encyclopedia,” 2017, https://en.wikipedia. 
org/w/index.php?title=Maglev&oldid=776151802. 

J. R. Reitz, “Forces on moving magnets due to eddy currents,” Journal of Applied Physics, 
41.5, 2067-2071, (1970). 

J. R. Reitz and L. C. Davis, “Force on a rectangular coil moving above a conducting slab,” 
Journal of Applied Physics, 43.4, 1547-1553, (1972). 

A. Singh, “Theory of eddy-current brakes with thick rotating disc,” Proceedings of the Institu- 
tion of Electrical Engineers, IET Digital Library, 124, no. 4, (1977) 

H. Hertz, Miscellaneous papers, (Macmillan and Co. Ltd, London, 1896), See Section 2, “On 
the induction in rotating spheres, 1880,” pp. 35-126. 

R. C. Youngquist et al., “A slowly rotating hollow sphere in a magnetic field: First steps to 
de-spin a space object.”, American Journal of Physics, 84, no. 3, 181-191, (2016) 

M. A. Nurge et al., “A thick-walled sphere rotating in a uniform magnetic field: The next step 
to de-spin a space object.” American Journal of Physics, 85, no. 8, 596-610, (2017) 

P. Moon and D. E. Spencer, Field Theory Handbook, 2nd Edition, (Springer-Verlag, Berlin, 
Heidelberg, 1988). 

P. Lorrain and D. Corson, Electromagnetic Fields and Waves, 2nd Edition, (W. H. Freeman 
and Co, New York, NY 1970), See Section 9..8: “Magnetic Field Calculation.” 

M. Schwartz, Principles of electrodynamics , (Dover Publications, Inc., Mineola, NY, 1987), pp. 
149-153. 


23 


15 


16 


17 


18 


19 


20 


21 


J. C. Maxwell, Electricity and Magnetism, (Clarendon Press, Oxford England. 1873), Vol. II, 
Secs. 694,695. 

A.H. Jaffey, “Solid angle subtended by a circular aperture at point and spread sources: formulas 
and some tables,” Review of Scientific Instruments, 25.4, 349-354, (1954). 

M. W. Garrett, “Solid angle subtended by a circular aperture,” Review of Scientific Instruments, 
25.12, 1208-1211, (1954). 

M. J. Prata, “Solid angle subtended by a cylindrical detector at a point source in terms of 
elliptic integrals,” Radiation Physics and Chemistry, 67.5, 599-603, (2003). 

J. C. Simpson et al., “Simple Analytic Expressions for the Magnetic Field of a Circular Current 
Loop,” NASA Technical Memorandum, NASA/TM-2013-217919, (2013). 

N. Praly et al., “Study on the eddy current damping of the spin dynamics of space debris from 
the Ariane launcher upper stages,” Acta Astro., 76, 145-153, (2012). 

F. Sugia et al., “Detumbling an uncontrolled satellite with contactless force by using an eddy 
current brake.” , 2013 IEEE/RSJ International Conference on Intelligent Robots and Systems, 
783-788, (2013) 

N.O. Gomez and S.J. I. Walker, “Eddy currents applied to de-tumbling of space debris: Analysis 


and validation of approximate proposed methods,” Acta Astro., 114, 34-53 (2015). 


24 


